A pan-genome of 69 Arabidopsis thaliana accessions reveals a conserved genome structure throughout the global species range.

Although originally primarily a system for functional biology, Arabidopsis thaliana has, owing to its broad geographical distribution and adaptation to diverse environments, developed into a powerful model in population genomics. Here we present chromosome-level genome assemblies of 69 accessions from a global species range. We found that genomic colinearity is very conserved, even among geographically and genetically distant accessions. Along chromosome arms, megabase-scale rearrangements are rare and typically present only in a single accession. This indicates that the karyotype is quasi-fixed and that rearrangements in chromosome arms are counter-selected. Centromeric regions display higher structural dynamics, and divergences in core centromeres account for most of the genome size variations. Pan-genome analyses uncovered 32,986 distinct gene families, 60% being present in all accessions and 40% appearing to be dispensable, including 18% private to a single accession, indicating unexplored genic diversity. These 69 new Arabidopsis thaliana genome assemblies will empower future genetic research.

Structural variation and DNA methylation shape the centromere-proximal meiotic crossover landscape in Arabidopsis.

BACKGROUND: Centromeres load kinetochore complexes onto chromosomes, which mediate spindle attachment and allow segregation during cell division. Although centromeres perform a conserved cellular function, their underlying DNA sequences are highly divergent within and between species. Despite variability in DNA sequence, centromeres are also universally suppressed for meiotic crossover recombination, across eukaryotes. However, the genetic and epigenetic factors responsible for suppression of centromeric crossovers remain to be completely defined. RESULTS: To explore the centromere-proximal meiotic recombination landscape, we map 14,397 crossovers against fully assembled Arabidopsis thaliana (A. thaliana) genomes. A. thaliana centromeres comprise megabase satellite repeat arrays that load nucleosomes containing the CENH3 histone variant. Each chromosome contains a structurally polymorphic region of ~3-4 megabases, which lack crossovers and include the satellite arrays. This polymorphic region is flanked by ~1-2 megabase low-recombination zones. These recombination-suppressed regions are enriched for Gypsy/Ty3 retrotransposons, and additionally contain expressed genes with high genetic diversity that initiate meiotic recombination, yet do not crossover. We map crossovers at high-resolution in proximity to CEN3, which resolves punctate centromere-proximal hotspots that overlap gene islands embedded in heterochromatin. Centromeres are densely DNA methylated and the recombination landscape is remodelled in DNA methylation mutants. We observe that the centromeric low-recombining zones decrease and increase crossovers in CG (met1) and non-CG (cmt3) mutants, respectively, whereas the core non-recombining zones remain suppressed. CONCLUSION: Our work relates the genetic and epigenetic organization of A. thaliana centromeres and flanking pericentromeric heterochromatin to the zones of crossover suppression that surround the CENH3-occupied satellite repeat arrays.

HEI10 is subject to phase separation and mediates RPA1a degradation during meiotic interference-sensitive crossover formation.

Reciprocal exchanges of DNA between homologous chromosomes during meiosis, or crossovers (COs), shuffle genetic information in gametes and progeny. In many eukaryotes, the majority of COs (class I COs) are sensitive to a phenomenon called interference, which influences the occurrence of closely spaced double COs. Class I COs depend on a group of factors called ZMM (Zip, Msh, Mer) proteins including HEI10 (Human Enhancer of Invasion-10). However, how these proteins are recruited to class I CO sites is unclear. Here, we show that HEI10 forms foci on chromatin via a liquid-liquid phase separation (LLPS) mechanism that relies on residue Ser70. A HEI10S70F allele results in LLPS failure and a defect in class I CO formation. We further used immunoprecipitation-mass spectrometry to identify RPA1a (Replication Protein A 1) as a HEI10 interacting protein. Surprisingly, we find that RPA1a also undergoes phase separation and its ubiquitination and degradation are directly regulated by HEI10. We also show that HEI10 is required for the condensation of other class I CO factors. Thus, our results provide mechanistic insight into how meiotic class I CO formation is controlled by HEI10 coupling LLPS and ubiquitination.

SCEP1 and SCEP2 are two new components of the synaptonemal complex central element.

The synaptonemal complex (SC) is a proteinaceous structure that forms between homologous chromosomes during meiosis prophase. The SC is widely conserved across species, but its structure and roles during meiotic recombination are still debated. While the SC central region is made up of transverse filaments and central element proteins in mammals and fungi, few central element proteins have been identified in other species. Here we report the identification of two coiled-coil proteins, SCEP1 and SCEP2, that form a complex and localize at the centre of the Arabidopsis thaliana SC. In scep1 and scep2 mutants, chromosomes are aligned but not synapsed (the ZYP1 transverse filament protein is not loaded), crossovers are increased compared with the wild type, interference is lost and heterochiasmy is strongly reduced. We thus report the identification of two plant SC central elements, and homologues of these are found in all major angiosperm clades.

Meiotic recombination is confirmed to be unusually high in the fission yeast Schizosaccharomyces pombe.

In most eukaryotes, meiotic crossovers (COs) are limited to 1-3 per chromosome, and are prevented from occurring close to one another by CO interference. The fission yeast Schizosaccharomyces pombe, an exception to these general rules, was reported to have the highest CO number per chromosome and no or weak interference. However, global CO frequency was indirectly estimated, calling for confirmation. Here, we used an innovative strategy to determine COs genome-wide in S. pombe. We confirmed weak CO interference, acting at physical distances compatible with the patterning of recombination precursors. We revealed a slight co-variation in CO number between chromosomes, suggesting that a limiting pro-CO factor varies between meiocytes. CO number per chromosome varies proportionally with chromosome size, with the three chromosomes having, on average, 15.9, 12.5, and 7.0 COs, respectively. This reinforces S. pombe's status as the eukaryote with the highest CO number per chromosome described to date.

HEIP1 is required for efficient meiotic crossover implementation and is conserved from plants to humans.

Crossovers (CO) shuffle genetic information and physically connect homologous chromosomal pairs, ensuring their balanced segregation during meiosis. COs arising from the major class I pathway require the activity of the well-conserved group of ZMM proteins, which, in conjunction with MLH1, facilitate the maturation of DNA recombination intermediates specifically into COs. The HEI10 Interacting Protein 1 (HEIP1) was identified in rice and proposed to be a new, plant-specific member of the ZMM group. Here, we establish and decipher the function of the Arabidopsis thaliana HEIP1 homolog in meiotic crossover formation and report its wide conservation in eukaryotes. We show that the loss of Arabidopsis HEIP1 elicits a marked reduction in meiotic COs and their redistribution toward chromosome ends. Epistasis analysis showed that AtHEIP1 acts specifically in the class I CO pathway. Further, we show that HEIP1 acts both prior to crossover designation, as the number of MLH1 foci is reduced in heip1, and at the maturation step of MLH1-marked sites into COs. Despite the HEIP1 protein being predicted to be primarily unstructured and very divergent at the sequence level, we identified homologs of HEIP1 in an extensive range of eukaryotes, including mammals.

High-frequency synthetic apomixis in hybrid rice.

Introducing asexual reproduction through seeds - apomixis - into crop species could revolutionize agriculture by allowing F1 hybrids with enhanced yield and stability to be clonally propagated. Engineering synthetic apomixis has proven feasible in inbred rice through the inactivation of three genes (MiMe), which results in the conversion of meiosis into mitosis in a line ectopically expressing the BABYBOOM1 (BBM1) parthenogenetic trigger in egg cells. However, only 10-30% of the seeds are clonal. Here, we show that synthetic apomixis can be achieved in an F1 hybrid of rice by inducing MiMe mutations and egg cell expression of BBM1 in a single step. We generate hybrid plants that produce more than 95% of clonal seeds across multiple generations. Clonal apomictic plants maintain the phenotype of the F1 hybrid along successive generations. Our results demonstrate that there is no barrier to almost fully penetrant synthetic apomixis in an important crop species, rendering it compatible with use in agriculture.

Joint control of meiotic crossover patterning by the synaptonemal complex and HEI10 dosage.

Meiotic crossovers are limited in number and are prevented from occurring close to each other by crossover interference. In many species, crossover number is subject to sexual dimorphism, and a lower crossover number is associated with shorter chromosome axes lengths. How this patterning is imposed remains poorly understood. Here, we show that overexpression of the Arabidopsis pro-crossover protein HEI10 increases crossovers but maintains some interference and sexual dimorphism. Disrupting the synaptonemal complex by mutating ZYP1 also leads to an increase in crossovers but, in contrast, abolishes interference and disrupts the link between chromosome axis length and crossovers. Crucially, combining HEI10 overexpression and zyp1 mutation leads to a massive and unprecedented increase in crossovers. These observations support and can be predicted by, a recently proposed model in which HEI10 diffusion along the synaptonemal complex drives a coarsening process leading to well-spaced crossover-promoting foci, providing a mechanism for crossover patterning.

The megabase-scale crossover landscape is largely independent of sequence divergence.

Meiotic recombination frequency varies along chromosomes and strongly correlates with sequence divergence. However, the causal relationship between recombination landscapes and polymorphisms is unclear. Here, we characterize the genome-wide recombination landscape in the quasi-absence of polymorphisms, using Arabidopsis thaliana homozygous inbred lines in which a few hundred genetic markers were introduced through mutagenesis. We find that megabase-scale recombination landscapes in inbred lines are strikingly similar to the recombination landscapes in hybrids, with the notable exception of heterozygous large rearrangements where recombination is prevented locally. In addition, the megabase-scale recombination landscape can be largely explained by chromatin features. Our results show that polymorphisms are not a major determinant of the shape of the megabase-scale recombination landscape but rather favour alternative models in which recombination and chromatin shape sequence divergence across the genome.

inGAP-family: Accurate Detection of Meiotic Recombination Loci and Causal Mutations by Filtering Out Artificial Variants due to Genome Complexities.

Accurately identifying DNA polymorphisms can bridge the gap between phenotypes and genotypes and is essential for molecular marker assisted genetic studies. Genome complexities, including large-scale structural variations, bring great challenges to bioinformatic analysis for obtaining high-confidence genomic variants, as sequence differences between non-allelic loci of two or more genomes can be misinterpreted as polymorphisms. It is important to correctly filter out artificial variants to avoid false genotyping or estimation of allele frequencies. Here, we present an efficient and effective framework, inGAP-family, to discover, filter, and visualize DNA polymorphisms and structural variants (SVs) from alignment of short reads. Applying this method to polymorphism detection on real datasets shows that elimination of artificial variants greatly facilitates the precise identification of meiotic recombination points as well as causal mutations in mutant genomes or quantitative trait loci. In addition, inGAP-family provides a user-friendly graphical interface for detecting polymorphisms and SVs, further evaluating predicted variants and identifying mutations related to genotypes. It is accessible at https://sourceforge.net/projects/ingap-family/.

Changing Gly311 to an acidic amino acid in the MATE family protein DTX6 enhances Arabidopsis resistance to the dihydropyridine herbicides.

In modern agriculture, frequent application of herbicides may induce the evolution of resistance in plants, but the mechanisms underlying herbicide resistance remain largely unexplored. Here, we report the characterization of rtp1 (resistant to paraquat 1), an Arabidopsis mutant showing strong resistance to the widely used herbicides paraquat and diquat. The rtp1 mutant is semi-dominant and carries a point mutation in the gene encoding the multidrug and toxic compound extrusion family protein DTX6, leading to the change of glycine to glutamic acid at residue 311 (G311E). The wild-type DTX6 with glycine 311 conferred weak paraquat and diquat resistance when overexpressed, while mutation of glycine 311 to a negatively charged amino acid (G311E or G311D) markedly increased the paraquat and diquat resistance of plants, whereas mutation to a positively charged amino acid (G311R or G311K) compromised the resistance, suggesting that the charge property of residue 311 of DTX6 is critical for the paraquat and diquat resistance of Arabidopsis plants. DTX6 is localized in the endomembrane trafficking system and may undergo the endosomal sorting to localize to the vacuole and plasma membrane. Treatment with the V-ATPase inhibitor ConA reduced the paraquat resistance of the rtp1 mutant. Paraquat release and uptake assays demonstrated that DTX6 is involved in both exocytosis and vacuolar sequestration of paraquat. DTX6 and DTX5 show functional redundancy as the dtx5 dtx6 double mutant but not the dtx6 single mutant plants were more sensitive to paraquat and diquat than the wild-type plants. Collectively, our work reveals a potential mechanism for the evolution of herbicide resistance in weeds and provides a promising gene for the manipulation of plant herbicide resistance.

The synaptonemal complex imposes crossover interference and heterochiasmy in Arabidopsis.

Meiotic crossovers (COs) have intriguing patterning properties, including CO interference, the tendency of COs to be well-spaced along chromosomes, and heterochiasmy, the marked difference in male and female CO rates. During meiosis, transverse filaments transiently associate the axes of homologous chromosomes, a process called synapsis that is essential for CO formation in many eukaryotes. Here, we describe the spatial organization of the transverse filaments in Arabidopsis (ZYP1) and show it to be evolutionary conserved. We show that in the absence of ZYP1 (zyp1azyp1b null mutants), chromosomes associate in pairs but do not synapse. Unexpectedly, in absence of ZYP1, CO formation is not prevented but increased. Furthermore, genome-wide analysis of recombination revealed that CO interference is abolished, with the frequent observation of close COs. In addition, heterochiasmy was erased, with identical CO rates in males and females. This shows that the tripartite synaptonemal complex is dispensable for CO formation and has a key role in regulating their number and distribution, imposing CO interference and heterochiasmy.

Modulation of evening complex activity enables north-to-south adaptation of soybean.

Soybean, a typical short-day crop, is sensitive to photoperiod, which is a major limiting factor defining its north-to-south cultivation range. The long-juvenile (LJ) trait is controlled primarily by the J locus which has been used for decades by soybean breeders to delay flowering and improve grain yield in tropical regions. The J gene encodes an ortholog of the Arabidopsis Evening Complex (EC) component EARLY FLOWERING 3 (ELF3). To identify modifiers of J, we conducted a forward genetic screen and isolated a mutant (eoj57) that in combination with j has longer flowering delay compared with j single mutant plants. Map-based cloning and genome re-sequencing identified eoj57 (designated as GmLUX2) as an ortholog of the Arabidopsis EC component LUX ARRHYTHMO (LUX). To validate that GmLUX2 is a modifier of J, we used trans-complementation and identified a natural variant allele with a similar phenotype. We also show that GmLUX2 physically interacts with GmELF3a/b and binds DNA, whereas the mutant and natural variant are attenuated in both activities. Transcriptome analysis shows that the GmLUX2-GmELF3a complex co-regulates the expression of several circadian clock-associated genes and directly represses E1 expression. These results provide mechanistic insight into how GmLUX2-GmELF3 controls flowering time via synergistic regulation of gene expression. These novel insights expand our understanding of the regulation of the EC complex, and facilitate the development of soybean varieties adapted for growth at lower latitudes.

The cohesin loader SCC2 contains a PHD finger that is required for meiosis in land plants.

Cohesin, a multisubunit protein complex, is required for holding sister chromatids together during mitosis and meiosis. The recruitment of cohesin by the sister chromatid cohesion 2/4 (SCC2/4) complex has been extensively studied in Saccharomyces cerevisiae mitosis, but its role in mitosis and meiosis remains poorly understood in multicellular organisms, because complete loss-of-function of either gene causes embryonic lethality. Here, we identified a weak allele of Atscc2 (Atscc2-5) that has only minor defects in vegetative development but exhibits a significant reduction in fertility. Cytological analyses of Atscc2-5 reveal multiple meiotic phenotypes including defects in chromosomal axis formation, meiosis-specific cohesin loading, homolog pairing and synapsis, and AtSPO11-1-dependent double strand break repair. Surprisingly, even though AtSCC2 interacts with AtSCC4 in vitro and in vivo, meiosis-specific knockdown of AtSCC4 expression does not cause any meiotic defect, suggesting that the SCC2-SCC4 complex has divergent roles in mitosis and meiosis. SCC2 homologs from land plants have a unique plant homeodomain (PHD) motif not found in other species. We show that the AtSCC2 PHD domain can bind to the N terminus of histones and is required for meiosis but not mitosis. Taken together, our results provide evidence that unlike SCC2 in other organisms, SCC2 requires a functional PHD domain during meiosis in land plants.

Mutation of DELAYED GREENING1 impairs chloroplast RNA editing at elevated ambient temperature in Arabidopsis.

Chloroplasts are important for plant growth and development. RNA editing in chloroplast converts cytidines (Cs) to uridines (Us) at specific transcript positions and provides a correction mechanism to restore conserved codons or creates start or stop codons. However, the underlined molecular mechanism is not yet fully understood. In the present study, we identified a thermo-sensitive mutant in leaf color 1 (tsl1) and found that TSL1 is allelic to DELAYED GREENING 1 (DG1). The missense mutation of DG1 in tsl1 mutant confers a high temperature sensitivity and impaired chloroplast development at an elevated ambient temperature in Arabidopsis. Subsequent analysis showed that chloroplast RNA editing at several sites including accD-1568, ndhD-2, and petL-5 is impaired in tsl1 mutant plants grown at an elevated temperature. DG1 interacts with MORF2 and other proteins such as DYW1 and DYW2 involved in chloroplast RNA editing. In vitro RNA electrophoretic mobility shift assay demonstrated that DG1 binds to RNA targets such as accD, ndhD, and petL. Thus, our results revealed that DG1 is important for maintaining chloroplast mRNA editing in Arabidopsis.

New perspectives on the plant PARP family: Arabidopsis PARP3 is inactive, and PARP1 exhibits predominant poly (ADP-ribose) polymerase activity in response to DNA damage.

BACKGROUND: Poly (ADP-ribosyl) ation (PARylation) is an important posttranslational modification that regulates DNA repair, gene transcription, stress responses and developmental processes in multicellular organisms. Poly (ADP-ribose) polymerase (PARP) catalyzes PARylation by consecutively adding ADP-ribose moieties from NAD+ to the amino acid receptor residues on target proteins. Arabidopsis has three canonical PARP members, and two of these members, AtPARP1 and AtPARP2, have been demonstrated to be bona fide poly (ADP-ribose) polymerases and to regulate DNA repair and stress response processes. However, it remains unknown whether AtPARP3, a member that is highly expressed in seeds, has similar biochemical activity to that of AtPARP1 and AtPARP2. Additionally, although both the phylogenetic relationships and structural similarities indicate that AtPARP1 and AtPARP2 correspond to animal PARP1 and PARP2, respectively, two previous studies have indicated that AtPARP2, and not AtPARP1, accounts for most of the PARP activity in Arabidopsis, which is contrary to the knowledge that PARP1 is the predominant PARP in animals. RESULTS: In this study, we obtained both in vitro and in vivo evidence demonstrating that AtPARP3 does not act as a typical PARP in Arabidopsis. Domain swapping and point mutation assays indicated that AtPARP3 has lost NAD+-binding capability and is inactive. In addition, our results showed that AtPARP1 was responsible for most of the PARP enzymatic activity in response to the DNA damage-inducing agents zeocin and methyl methanesulfonate (MMS) and was more rapidly activated than AtPARP2, which supports that AtPARP1 remains the predominant PARP member in Arabidopsis. AtPARP1 might first become activated by binding to damaged sites, and AtPARP2 is then poly (ADP-ribosyl) ated by AtPARP1 in vivo. CONCLUSIONS: Collectively, our biochemical and genetic analysis results strongly support the notion that AtPARP3 has lost poly (ADP-ribose) polymerase activity in plants and performs different functions from those of AtPARP1 and AtPARP2. AtPARP1, instead of AtPARP2, plays the predominant role in PAR synthesis in both seeds and seedlings. These data bring new insights into our understanding of the physiological functions of plant PARP family members.

Amino acid based de Bruijn graph algorithm for identifying complete coding genes from metagenomic and metatranscriptomic short reads.

Metagenomic studies, greatly promoted by the fast development of next-generation sequencing (NGS) technologies, uncover complex structures of microbial communities and their interactions with environment. As the majority of microbes lack information of genome sequences, it is essential to assemble prokaryotic genomes ab initio aiming to retrieve complete coding genes from various metabolic pathways. The complex nature of microbial composition and the burden of handling a vast amount of metagenomic data, bring great challenges to the development of effective and efficient bioinformatic tools. Here we present a protein assembler (MetaPA), based on de Bruijn graph searching on oligopeptide spaces and can be applied on both metagenomic and metatranscriptomic sequencing data. When public homologous protein sequences are involved to guide the assembling procedures, MetaPA assembles 85% of total proteins in complete sequences with high precision of 83% on real high-throughput sequencing datasets. Application of MetaPA on metatranscriptomic data successfully identifies the majority of actively transcribed genes validated in related studies. The results suggest that MetaPA has a good potential in both metagenomic and metatranscriptomic studies to characterize the composition and abundance of microbiota.